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SHARP ERROR BOUNDS FOR PIECEWISE POLYNOMIAL 
APPROXIMATION: REVISIT AND APPLICATION TO ELLIPTIC 
PDE EIGENVALUE COMPUTATION 

HEHU XIE* AND ZHIMIN ZHANGt 

Abstract. In this paper, we revisit approximation properties of piecewise polynomial spaces, 
which contain more than Vr—i but not Vr- We develop more accurate upper and lower error bounds 
that are sharper than those used in literature. These new error bounds, especially the lower bounds 
are particular useful to finite element methods. As an important application, we establish sharp 
lower bounds of the discretization error for Laplace and 2m-th order elliptic eigenvalue problems 
in various finite element spaces under shape regular triangulations, and investigate the asymptotic 
convergence behavior for large numerical eigenvalue approximations. 
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1. Introduction. Standard error estimates for numerical methods are conven¬ 
tionally presented as upper bounds. For example, the error estimate for the finite 
element method is usually written in the following way 

\\U- Uh\\l <Ch'^~^\\u\\r, (1.1) 

where Uh denotes the finite element approximation (of polynomial degree less than 
r) to the exact solution u. However, this error estimate is not so delicate and sharp 
for some finite elements such as the tensor product and intermediate families. The 
first aim of this paper is to establish some more delicate and sharp upper bounds for 
piecewise polynomial approximation errors on tensor product meshes by giving a new 
norm equivalence theorem. 

The analysis of lower bound error estimates have been studied in the literature 
[SEiiniiisiiiiiiii]. However, none of them are particularly for eigenvalue approxi¬ 
mation. In particular, when eigenvalue problems are involved, these upper bounds are 
tied with the magnitudes of the underlying eigenvalues. For the same finite element 
method, larger eigenvalues always have larger error bounds, the fact was described 
in Strang-Fix’s 1973 classical book [16]. Inspired by some recent studies on effi¬ 
ciency of finite element approximation for different eigenvalues [20] . this work is also 
devoted to establishing some lower bound error estimates for finite element approxi¬ 
mation of eigenvalue problems. Then, the second aim of this paper is to derive lower 
error bounds for eigenvalue approximations by finite element methods with contin¬ 
uous or discontinuous piecewise polynomial functions under reasonable assumptions 
on regularities of eigenfunctions and triangulations. As an application, we show the 
asymptotic convergence behavior for the large numerical eigenvalue approximations. 
This knowledge can help us to estimate the number of reliable numerical eigenvalues 
that convergence at a predetermined rate, and thereby provide a solid theoretical 
foundation for the issue raised and discussed in EQ]. 
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Towards this end, we revisited piecewise polynomial approximations and obtain 
more delicate error estimates in terms of higher-order derivatives. Then combining 
norm estimates of an eigenfunction in terms of the underlying eigenvalue leads to a 
lower bound of the approximation error for the eigenvalue problem in the following 
form (see Sections 3 and 4): 

\\u-Uh\\j,p,h>CX^h‘^~^, 0<j<r, 

where Uh is the finite element approximation (of polynomial degree less than r) to the 
Laplace eigenvalue problem, and C is a constant independent of the mesh size h and 
the eigenvalue A. 

The outline of the rest of the paper goes as follows. Section 2 is devoted to 
introducing the notation and a new norm equivalence theorem. Then some upper 
bounds of the error by piecewise polynomial approximation are given in Section 3. 
Section 4 is for lower bounds of the discretization error of the Laplace eigenvalue 
problem by finite element methods. The lower bounds of the discretization error of the 
2TO-th order elliptic eigenvalue problem are presented in Section 5. Some concluding 
remarks are given in the last section. 

2. A norm equivalent theorem: Revisit. In this section, we first introduce 
some notations, and then state a new norm equivalence theorem, from which we derive 
new upper bounds of the error by piecewise polynomial approximation. 

Here we assume that fl C 7^" (n > 1) is a bounded polytopic domain with 
Lipschitz continuous boundary dfl. Throughout this paper, we use the standard 
notation for the usual Sobolev spaces and the corresponding norms, semi-norms, and 
inner products as in [mms]. 

n 

A multi-index a is an n-tuple of non-negative integers with length |q;| = 

i—1 

The derivative is then defined by 


/ i9 N“1 

• ( ^ ) 


\dXn^ 


For a subdomain G of fl, the usual Sobolev spaces W"^’P{G) with norm || • \\m,p,G and 
semi-norm | ■ \m,p,G are used. In the case p = 2, we have H'^{G) = VF"‘’^(G') and the 
index p will be omitted. The L^-inner product on G is denoted by (•, ■)a. For G C H 
we write G CC O to indicate that dist(9fl, G) > 0 and meas(G) > 0. 

We introduce a face-to-face partition 7^ of the computational domain H into 
elements K (triangles, rectangles, tetrahedrons, bricks, etc.) such that 

y A 

Ken 

and let £h denote a set of all (n — l)-dimensional facets of all elements K G Th- Here 
h := maxxeTh and Hk = diam K denote the global and local mesh size, respec¬ 
tively [Hi. We also define = {K GTh and K C G} and he = max^^g^-G Kk- A 
family of partitions Th is said to be shape regular if it satisfies the following condition 

0 : 


3cr > 0 such that hxlrx > cr MK G Th, 
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where tk is maximum diameter of the inscribed ball in K G Th- A shape regular 
family of partitions Th is called quasi-uniform if it satisfies [3 [5] 

3/? > 0 such that ma.x{h/hx, K G Th} < /3. 


A finite element space Vh of piecewise polynomial functions is constructed on Th as 

Vh = {v.v\KGV^,,A, yKGTh}, (2.1) 

where T^used denotes the polynomial space associated with multi-index set Indused 

■Pused = I'C : P = ^ 

aelnduaed 


We define multi-index sets 

Indr = {a : la] < r}, Indr,rest = Indr\Indused, 

and assume that Indr-i C Indused and Indr ^ Indused- We further define the following 
broken semi-norm for v € U 14 with G C Q 


and 


1 


\v\j,p,G,h ■— 1 

s L 

X \D^v\^dK\ , 

1 < p < oo 



4l=i / 



klj.oo,G,h := niax |u|j,oo.A. 
Then the corresponding norm can be defined by 


j,p,G,h ■— [ ^ 


\i^p^G^h j ’ 


U=0 


and 


•— max 

0<i<j 


We will drop G when G = fl. Throughout this paper, the symbol C or c (with or 
without subscript) stands for a positive generic constant which may attain different 
values at its different occurrences and which is independent of the mesh size h. 

Now, we state a more general norm equivalence theorem, of which the conventional 
norm equivalence theorem is a special case. 

Theorem 2.1. Let LI C 7?." be a Lipschitz domain. For j = 1, 2,..., J, Assume 
that fj : W’’+^’^(n) —>■ TZ with r > 1 and 1 < p < oo are semi norms on 
satisfying two conditions: 

(HI) 0 < f,{v) < C||u||,+i,p.n, Vu G W^^P{n), l<j<J. 

(H2) If V G Vused 0 ‘nd fj{v) = 0 for 1 < j < J, then u = 0. 

Then the quantity 


J 

Ikll = XI + + X P'^^llo.p.n 

yelndr.reat i=l |7|=r-|-l 


( 2 . 2 ) 
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or 

lkll=f E P"Hlg.p.n + E/.(d^+ E ll^"Hlo.p.n) (2.3) 

y7Glndr,rest i=l |7|=r+l J 

defines a norm on which is equivalent to the norm ||i;||j.+i_p_Q. 

Proof. We prove that the quantity (12.21) defines a norm on equivalent 

to the norm ||?^||r+i,p,a. The statement on the quantity (12.31) can be proved similarly 
or by noting the equivalence between the two quantities (12.21) and (12.31) . 

By the condition {HI), we know that for some constant C > 0, the following 
estimate holds 


Ikll < C||^;||,+i,p,a V« e W^+^^P{n). (2.4) 

So we only need to prove that there exist another constant C > 0 such that 

Ikllr+I.p.a < c||u|| Vu e (2.5) 

We argue by contradiction. Suppose the inequality is not true. Then we can find a 
sequence {ufe} C W’'+^’^'(n) with the properties 

lkfe||r+i.p,a = 1, (2.6) 

Ikfell < I (2.7) 

for fc = 1, 2,.... From (|2.7I) . we can see that as fc —?► oo, 

E (2.8) 


and 


E \\D^^kh,p,a ^ 0 

|7l=»'+l 


(2.9) 


fj{vk)^0, 1<3<J. ( 2 . 10 ) 

Since {vk} is a bounded sequence on W''^^’^{^1) from the property (12.61) and W’'+^’P(fl) 
can be embedded compactly into VF’’’P(n), there is a subsequence of the sequence {wfc}, 
still denoted as {vk}, and a function v G W^’^{V.) such that 

Vk ^ V in W^’^{D,), as k ^ oo. (2-11) 

This property and dH-dlsi), together with the uniqueness of a limit, imply that 

Vk ^ V in the norm ||u||r+i,p.n, 


E llT^^^llo.p.n + E 

7elndr,reat |7|=r + l 


as k ^ oo, and 
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^ I X] ll^^^fcllo,p.n + X Il-O^^fello,p,n 1 = 0- 

y7elndr,reat |7|=’'+1 / 

We then conclude that v € 7^used(^^)- On the other hand, from continuity of function- 
(12.7L we find that 

fj{v) = lim fj{vk) = 0. 

k—^oo 

Using the condition { H 2 ), we see that n = 0, which contradicts the relation that 

Ikllr+I.p.n = lim ||nfc||r+i,p,n = 1- 
k—¥oo 

Hence the inequality (12.51) holds which means the norm || • || is equivalent to the one 
II • ||r+i.p.n and the proof is completed. □ 

In the error analysis for the finite element method, we also need an inequality 
involving the norm of the Sobolev quotient space 

H = W^+^’P{n)/rr,used = {M : H = {^; + q\q e rr,used} j(?.12) 

where Vr,used denotes the polynomial space associated with multi-index set Indr^used = 
Indj. n Indused. 

Any element [n] of the space V is an equivalence class, the difference between any 
two elements in the equivalence class being a polynomial in the space 'Pr.used- Any 
[u] is called a representative element [u]. The quotient norm in the space V is defined 
to be 


IIMIk= inf ||?; + g||^+i,p.n. 

<?ePr-,used(i^) 

Corollary 2.2. Let H C he a Lipschitz domain. Then the quantity 

ll^^'*^llo.p.a +E| 7 |=r-n P^i^llo.p.a (1 < P < oo) is a semi-norm on V, 
equivalent to the quotient norm ||[n]||y. It means 

Cl I X l|C’'^^llo,p.n + X P'^^llo.p.o) < inf ||n-hg||r+i,p,n 
V76ind.,„. |7|=.+i J 

< 02 ! Y. \\D^^\\o,p,n+ Y P^^llo.P.o), (2.13) 

y76lndr,rest |7|=’'-I-1 / 

where Ci and C 2 are two constants independent of v. 

Proof. Obviously, for any [n] G V and any v € [u], we have 

IIHIk= inf ||ii + 9l|r-m,p.n > Y P^^llo.p.n + Y l|C’'^^llo.P.n- 

96W,used(f2) d I I , , 

7elndr,rest laHr-l-l 

Thus we only need to prove that there is a constant C, depending only on H, such 
that 


inf ||ii + 9||r-i-i,p.a < C E l|P''dlo,p,a + E l|P’dl»,p.ii b.l4) 
9eP.,u,ed(n) V7Glnd.,„.t |7l=--+l 
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Denote N = dim(7^j. used(^))- Define N independent linear continuous functionals on 
'Pr,used: the Continuity being with respect to the norm of Now we can use 

Theorem [Q to obtain the following inequality 


||r+l,p,0 C I i: +1: iic-'.iio,p,ii 

v7Glndr,peat | 7 |='r+l 


E 


D'^vdD, 


, Vu e (2.15) 


Replacing v hy v + q and noting that q = 0 for 7 € Ind^^est U Ind^+i, we have 


Ik + 9llr-+i,p,n < C 


E ll^'^klo.p.a 

7elndr,rest 


|7l=»’+l 



Vu e lR’'+hP(0) and Vg e P^.used- 


Now we construct a polynomial q G Pj-.used satisfying 


(2.16) 


/ D'^{v + q)dQ = 0 , for 7 G Indr.used- (2.17) 

Jo, 

This can always be done: set 7 £ Indr.used, then D'^q equals 7 ! = 71 ! ■ • ■ 7 ^! times the 
coefficient of := so the coefficient can be computed by using (12.1711 . 

Having found all the coefficients for terms in T^r-.used, we set 7 G Indr,used, and use 
(12.1711 to compute all the coefficients for terms of Vr,used- Proceeding in this way, we 
obtain the polynomial q satisfying the condition (I2.17|l for the given function v. 

With g = g in (12.1611 . we have 


inf ||u + g||r-+i,p,n < lk + g||r+i,p,a 

9e'Pr,used(f2) 

<cl P'^klo.p.o+ E ll^^klo.p.a 

y7elndr,rest 171=^+! 

This is the desired result (12.1411 and the proof is complete. □ 

Corollary 12.21 is more delicate than the standard norm equivalence theorem. In 
some special meshes, this new estimate can lead to new and sharper error estimates 
for the piecewise polynomial approximation. 

3. New upper error bounds for piecewise polynomial approximation. 

In this section, as applications of the new norm equivalence theorem stated in Theorem 
[231 we give some error estimates for the tensor product and intermediate finite element 
spaces. 

Theorem 3.1. Assume u G W’'“'"^’^(f2) and the mesh is tensor-product type with 
the mesh size hK.i,''' j in each direction and the edges paralleling to the axes 
for each K G Th- IPe use the finite element space 14 defined in a to approximate 
the function u. Then the following upper bound of the error holds 



inf \\D^{u 

Vh&Vh 


'>^h)\\o,p,Q,h < C 


4 E 

K&Th 7Glndr,rest 
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E E 

K^Th |7|=r+l 


P(7-a) 

K 




0,p,-fC 


, (3.1) 


where . 

Proof. The proof can be given by combing Corollary 12.21 Indr,used Q Indused and 
the standard scaling argument (cf. [HIS]). □ 

Theorem O can be used to derive sharper (than conventional) error estimates 
for tensor product and intermediate finite element spaces [iiziii]. The tensor product 
finite element space on tensor product meshes is defined as follows 0[2l[8]: 


Vh = : VklK G Qr-l, VitT e 7 )i|. 


(3.2) 


Here and hereafter Qk denotes the space of polynomials with a degree no more than k 
for each variable, i.e. Qk = {x°‘ with ai < k for i = 1, • • • ,n}. For example, when Vh 
is the biquadratic element with r = 3 and n = 2, we know Inds.rest = {(3,0), (0,3)} 
and we have the following error estimates 


inf ||u U?i||o<C( f^%,l\\9xxxU\\Q + h^j^2\\9yyy'U'\ 


Vh&Vh 


\KeTh 


+C'I E E 

^KgTh | 7|=4 


(3.3) 


inf ||a:,(u-u,,)||o < C h%^\\dxxxu\\o + hYh%,2\\9yvvu\\ 


+C I E E 


0,2 


iKGTh |7N4 


(3.4) 


inf ||5y(u-U,,)||o < C Y >T'K,lhj.^W\dx:xxU\\l + h% 2 \\dyyyU\\ 


I E E V 

^KGTh | 7|=4 


2 ( 7 -( 0 . 1 ))||^ 7„||2 


l0,2 


(3.5) 


which are different from the traditional error estimates. 

The intermediate family of the second type on rectangular meshes is defined as 
follows Ii[2T]: 


Vh = ^Vh ■ V/iIk ^ 'Pr+l n Qr-l, ViF € 7)i |, 


(3.6) 


where Vk denotes the space of polynomials with degree no more than fc, i.e. Vk = {a;“ 
with |a| < k}. For example, when r = 4, we have Ind 4 ,rest = {(4,0), (0,4)} and the 
following error estimates 


inf ||u U?,||o < C [ 'Y ^%,i\\9xxxxU\\q + h^2\\9yyyyU\\Q 

\iYr, / 
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I E E 

.Ken |7l=5 


(3.7) 


inf ||i9a;(u 'f/i)||o ^ C" ( X! i^^, 2 \\^yyyy'^\ 

\£tn 


+c\ E , 

.Ken |7l=5 


(3.8) 


ll^y(^ ^?i)||o ^ C" ( X! ^%.l^K 2\\^^^^^'^\\o ^%.2\\9yyyy'^\ 

\Ken 


+GI E E 

, i^eT;. |7|=5 


2(7-(0.1)) 


\\D-^u\\l, 


(3.9) 


Theorem ED can also be used to give a new upper bound error estimate for the 
serendipity finite element families Enaiz]: 


Vfi = : Vh\K s Sr, \fK G 7)i|, 


(3.10) 


where Sr = Pr-i+spanja;’' ^y, xy^ ^}. For example, when r = 3, we know Inda^rest = 
{(3,0), (0,3)} and the following error estimates hold 


inf ||u - u^^llo < C' 'Y' h%j^\\dxxxu\\o 
\i^n ’ 




'yyy'^Wo 


+C' I E E 

\Ken |7|=4 


2 


(3.11) 


inf ||a,,(M-nft)||o < C h%^\\dxxxu\\l + hYihK,2\\9yyyu\\ 

\i^n 


+c I E E 


0,2 5 


kK^Th |7N4 


(3.12) 


iiif i^?i)||o < C* [ X! ^k,i^k,2\\^^^^'^\\q ^K,2\\^yyy'^\\ 

\i^n 


+c I E E 


0,2 5 


^Ken |7l=4 


(3.13) 


which are different from the standard error estimates. 

Compared with the traditional upper bound (ini), the above new bounds are 
sharper in some sense. Theorem 13.11 can also been used to deduce some more delicate 
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and sharper error estimates for other finite element families. To end this section, we 
would like to point out that this type of error estimates (with some difference) have 
also been given and proved in [71 Chapter 4.5] by the average Taylor polynomials tech¬ 
nique. Nevertheless, we need more delicate theory in our analysis and the approach 
based on the quotient space argument in this paper serves the purpose and requires 
minimum regularity assumption. 

4. Lo-wer bounds for finite element approximation of Laplace eigen¬ 
value problem. Another important issue in this paper is to derive the lower bounds 
of the discretization error for the Laplace eigenvalue problem by the finite element 
method. The special feature here is that the discretization error will be expressed in 
the power function of the corresponding eigenvalue. Combined with the upper bounds 
of the discretization error Him, we are able to estimate the number of reliable nu¬ 
merical eigenvalues at the predetermined convergence rate. For more information and 
examples, please refer to |20j . Here we provide the corresponding theoretical analysis. 
Here we are concerned with the following Laplace eigenvalue problem: 

Find (A,m) such that ||u||o = 1 and 

f —Alt = Alt in H, 

I „ = 0 o„aa 

Based on the partition l~h on f2, we can define a suitable finite element space Vu 
(conforming or nonconforming for the second order elliptic problem) as in (12.11) with 
piecewise polynomials of degree less than r. We solve the eigenvalue problem (14.1|) 
by the finite element method. The finite element approximation {Xh,Uh) GTZxVh of 
(03 can be defined as follows: 

Find {Xh,Uh) €TZxVh such that ||?tti||o = 1 and 

ah{uh,Vh) = Xhiuh,Vh) VvhGVh. (4.2) 

First, it is well known that the following upper bound of the discretization error (see 
[HIHlIle]) holds 

\\u-Uh\\e,p,h<Ch‘^-^X^^-^'>/^, 0<£<1, 0 < s < r, (4.3) 

where 1 < p < oo and the constant C is independent of the mesh size h and the 
eigenvalue A. In the rest of this section, we will prove is also a lower 

bound of discretization error ||u — Uh\\t,p,h under some suitable conditions. 

Lemma 4.1. For the eigenvalue problem eip. the following property holds 

||Vl'51-L§lALiJu||o,p,G = AL5J||vl'5’l“L51y||g^p_g,^ ( 44 ) 

where G CC |"r] denotes the smallest integer not less than r and [rj is the biggest 
integer not larger than r. 

Proof. From the eigenvalue problem (14.11) . we have the following equation 
V^il-Lil ALiJu = aLSJ 

which leads to the desired result (14.41) and the proof is complete. □ 

From Lemma im we state the following lower bound results of the discretization 
error. 

Theorem 4.2. Assume that = 0 for all K G Th, Vh G Vh, where 

the partition Th is quasi-uniform and shape regular. If there is a domain G CC H 
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such that ||w||o,p,G > 0; then for any given exact eigenpair (A,u); its corresponding 
eigenpair approximation (A/j, Uh) €TZxVh in satisfies the following lower bound 

of the discretization error 

\\u-Uh\\j,p,h>Ch'^~^\^, 0<j<r, (4.5) 

where the mesh size h is small enough, 2 < p < oo and C is a positive constant 
independent of X and the mesh size h. 

Furthermore, if the family {Vh} of partitions is only shape regular, then for any 
given exact eigenpair {X, u), its corresponding approximation (Xh,Uh) has the following 
lower bounds of the discretization error 




>CX^, 0<j<r 


(4.6) 


and 


h-K Ir I 0<j<r, 

yK&T^ 


(4.7) 


where the mesh size h is also small enough, 2 < p < oo, p < q < oo, C are positive 
constants independent of X and the mesh size h. 

Proof. First, it is easy to see that the eigenfunctions of problem 63 cannot be 
polynomials of bounded degree on any subdomain G CC fl. It means that we also 
have ||Vm||o,p,g > 0. 

Now, we present the derivation for the two cases of the positive integer r. Toward 
this end, let 11)) denote a suitable piecewise Vr interpolation operator (e.g., piecewise 
L^-projection). 

In the first case, r is even. Then we have 

||A2u||o,p.G = ll^^(u — Vh)\\o,p,G,h 

< ||A5(u - Ulu)\\o,p,G,h + ||Ai(n> - Vh)\\o,p,G,h 

< Ch^\\u\\r+S,p,G + Ch^-^\\nlu - Vh\\j,p,G.h 

< Ch^\\u\\r+S,p,G + Ch^-^\\nlu - u\\j,p,G,h 

+Ch^~^\\u - Vh\\j,p,G,h 

< Ch^\\u\\r+S,p,G + Ch^~'~\\u — Vh\\j,p,G,h- 


It means 


Ch^-^\\u - Vh\\j,p,G > C\\Aiu\\o,p,G - Ch^\\u\\r+S,p,G- 
Together with Lemma l4.11 when h is small enough, we have 

Ch^-''\\u-Vhhp,G,h > G||A5n||o,p.G = GA5||n||o.p.G > CX^ . 
In the second case, the integer r is odd. Then we have 
II VA“2-7i||g = II VAt“ [u — f?i)||o.p.G 


(4.8) 


(4.9) 
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< IIVA^ {u - n»||o.p,G + IIVA^ (n> - z;^)||o.p.G 

< Ch^\\u\\r+S,p,G + Ch^-^\\Ulu - VhWj,P,G 

< Ch^\\u\\r+S,p,G + Ch^-^wnlu - u\\j,p,G 

+ Ch^-^\\u-Vh\\j,p,G 

< Ch^\\u\\r+S,p,G + Ch^-^\\U - VhWj,p,G 

which means 

Ch^~^\\u — Vh\\j,p,G > C'l|VA“ 2 -y||p_p g _ C'/l^||'u||r+5.p,G- 

Similarly, together with Lemma 14.11 when h is small enough, the following estimate 
holds 

Ch^-^\\u-Vh\\j,p,G > C||VA^u||o.p,g = CA^IIVmIIo.p.g 

>qo|A^||Vn||o.G>C'A5, (4.10) 

where we used ||Vu||o,n = A^^^ and the constant C depends on || Vu||o,g/II Vu||o,a. 
Combining (14.91) . (14.101) and the arbitrariness of Vh leads to 



Vh&Vh 


'^h\\j,p,G,h 


> CA^. 


Together with the following property 

IIm ~ Uh\\j,p,h ^ ||m ~ Uh\\j,p,G,h ^ iq ~ Vh\\j,p,G,h 

h^-3 ~ - vhdVh h^-3 


we obtain the desired result (H3]) . The error bounds (ITO and dHH) can be proved 
similarly. □ 

Remark 4.1. The interior regularity result u € for a subdomain 

G CC n and 5 > 0 for the Laplace eigenvalue problem o can be obtained from m 
Theorem 8.10 and P.214 with some recursive discussion] when the right-hand side f 
is sufficiently smooth. 

Now, we present some conforming and nonconforming elements whose discretiza¬ 
tion errors attain lower bounds in the forms of (14.51) , (14.61) , and (14.71) with the help of 
Theorem 14.21 

First we can obtain the lower bound results for the standard Lagrange type ele¬ 
ments 

Vh = [vh\K e VeiK) or Qi(i^), ViF e %]■ (4.11) 

From Theorem the lower bound results in this section hold with r = i 1 for 
Vt{K) case and r = 2 for Qi{K) case. 

Then it is also easy to verify the lower bound results for nonconforming elements 
Crouzeix-Raviart (CR) and Qi rotation (Q™*): 

• The CR element space, proposed by Crouzeix and Raviart [9], is defined on 
simplicial partitions by 

Vh = [v&L^{n)-. v\k(1Vi{K), 

J v\K^ds = J v\K 2 ds if RTi n R :2 = G 


The lower bound result holds with r = 2 and 7 S Ind 2 . 
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• The Q™* element space, proposed by Rannacher and Turek and Arbogast 
and Chen [5], is defined on n-dimensional block partitions by 


14 = ji’ e L^{n) : v\k £ QRot(A'), 



where Q ^ ot ^ K ) = Vi { K )+ span { x ^ — \ l<*<n —l}. The lower bound 

result holds with r = 2 and 7 with 7 j = 1 , 7 ^- = 1 , 1 < i < j < n. 

Comparing (H3)) and g3D . we see that all lower and upper bounds of the above 
examples are sharp if the solution is smooth enough such that the upper bounds hold. 

Now, we investigate the number of reliable eigenvalue approximations obtained 
by the finite element method via a rigorous proof for the procedure outlined in | 20 | . 
The techniques for the analysis are the lower bound result (14.51) and the asymptotic 
behavior for the j-th eigenvalue of the Laplace operator. In the following analysis, we 
need the following estimate which comes from m Chapter 6 , Section 6.3]: 


llMllr+ 1,2 < C'll^ ^ ^llo,2 < CX 2 . 


(4.12) 


Lemma 4.3. Assume that the mesh Th is quasi-uniform and we use the conform¬ 


ing finite element method. Under the condition of Theorem \4.^ we have the following 
lower bound for the numerical eigenvalue relative error 




> Ch^-^\\A^u\\o.2,G - Ch’'X^ 

> Ch^-^Xi - Ch^X"^, yvh G Vh. 


(4.14) 


The desired result (14.131) can be obtained by combing (14. 5|) and the following result 



When r is odd, we can also deduce the desired result (14.131) in the similar way. □ 


In [T^, Weyl gives the following asymptotic behavior for the j-th eigenvalue of 
the Laplace operator: 


2 



(4.15) 


where a;„ = /T{1 + n/2) denotes the volume of the unit ball in . 

Based on the above preparation, we come to prove the results in [20j . 

Theorem 4.4. Let N denote the total degrees of freedom in the finite element 


space Vh- Assume the conditions of Lemma\4.3\ are satisfied andjg = for 9 G [0,1). 


Then the relative error of the jg-th numerical eigenvalue has the following lower bound 



(4.16) 
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and the absolute error of the jg-th numerical eigenvalue has the following lower bound 
K,h - K > C{r - lf 0 rh^r(l-e)- 2 ^ 


when N is large enough. 

On the contrary, the number jg of the “reliable” numerical eigenvalues with rel¬ 
ative error of converging at the rate for 9 S (0,1], has the following 

estimate 


jg < C{r - (4.18) 

Furthermore, the number jg of the “reliable” numerical eigenvalues with absolute error 
of Xjg converging at the rate g ^ (q^ following estimate 

jg <C{r— ?-. (4.19) 

Proof. First, the piecewise polynomial space of degree r — 1 has the total degree 
of freedom N = 0{{r — I)"/*””) which means 

h = 0{N-i{r-l)-^). (4.20) 

Then the desired results (|4.16|) and (I4.17p can be derived by combining (|4.13|) . (|4.15p . 
(I4.20p and the following estimate 

jg < C{r — l)“^iV"^®“^^ ^ 1, when 9 < I and N is large enough. (4.21) 

Bv usine the lower bound results (|4T3|) and (|4.15p and along the wav in m , we 
can give the proof for (|4.18p . From (j4.13p and (|4.15p . we have 

Then the desired convergence order leads to the following estimate 

h«(2('-i)) >C'(j^h2)’-'(l-jth2). (4.22) 

Combining (14.2011 . (14.2111 and (14.2211 leads to the desired number jg satisfying the 
inequality (I4.18P and (14.1911 can be proved similarly. □ 

5. Lotver bounds for finite element approximation of 2m-th order el¬ 
liptic eigenvalue problem. We consider similar lower bounds of the discretization 
error for the 2m-th order elliptic eigenvalue problem by the finite element method. 
This is a natural generalization of the results in Section 0] 

The 2m-th order Dirichlet elliptic eigenvalue problem for a given integer to > 1 
is defined as 


(_l)"^A™u 

u 

l|w||o 


Xu in ft, 

0 on dfl and 0 < j < to — 1, 

1 , 


(5.1) 


where v denotes the unit outer normal. The corresponding weak form of problem 
(EB) is defined as follows: 
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Find (A,m) GTZx such that ||m||o = 1 and 

a{u^v) = X{u,v) Vv £ (5-2) 

where 


a{u,v)= [ y^uS/^^vdn. 

Jn 

Based on the partition 7^ of 12, we build a suitable finite element space Vh (con¬ 
forming or nonconforming for the 2m-th order elliptic problem) with a piecewise 
polynomial of a degree less than r and define the corresponding discrete eigenvalue 
problem in the finite element space: 

Find {Xh,Uh) GTZxVh such that ||w?i||o = 1 and 

ah{uh,Vh) = Xh{uh,Vh) Vvh£Vh, (5.3) 


where 

ah{uh,Vh)= ^ j 

Ken 

Now we conduct the corresponding lower-bound analysis for the eigenvalue prob¬ 
lem H5.2I) . Similarly to Lemma 14.11 the following estimate for the eigenfunction norm 
by the eigenvalue holds. 

Lemma 5.1. For the eigenvalue problem the following property holds 

||V™A™^u||o,p.g = A^I|V™u||o.p,g, (5.4) 

where G CC 12, i = 0,1 and £ = 0,1, • • •. 

From Lemma l5.ll we state the following lower bound results of the discretization 
error. 

Theorem 5.2. Let r = mi + 2m£. Assume that = 0 for i = 0,1 

and i = 0,1, ■ ■ ■ for all K £ Th, Vh £ Vh, where the partition Th is quasi-uniform and 
shape regular. If there is a domain G CC 12 such that ||m||o,p,g > 0, then for any given 
exact eigenpair (X,u), its corresponding eigenpair approximation {Xh,Uh) £ TZ x Vh 
in JOI) satisfies the following lower bound of the discretization error 

\\u-Uh\\j,p,h>Gh'^~^X^, 0<j<r, (5.5) 


where the mesh size h is small enough, 2 < p < oo and C is a positive constant 
independent of u, X and the mesh size h. 

Furthermore, if the family {Tf} of partitions is only shape regular, then for any 
given exact eigenpair {X, u), its corresponding approximation {Xh,Uh) has the following 
lower bounds of the discretization error 


E ^ 


LPU-r) 

'■iC 


\U - Uh 



> GX-- 


0<J< 


and 


,p{U-i")+n( 

E 


D) 


U-Uh 



(5.6) 


> GX^ 


0 < j < r. 


(5.7) 
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where the mesh size h is small enough, 2 < p < oo, p < q < oo, C are positive 
constants independent of u, X and the mesh size h. 

Proof. First, it is also easy to see that the eigenfunctions of problem (15.11) cannot 
be polynomials of bounded degree on any subdomain C? CC fl. It means that we also 
have ||V™u||o,p,G > 0. 

Now, we present the derivation for the two cases of the integer i. In the first case, 
i = 0 and r = 2mi. Then we have 

l|A'"^u||o.p.G = l|A™^(u-u,)||o,p.G.h 

< ||A™^(u - n»||o.p.GA + ||A"*'(n> - Vh)\\o,p,G,h 

< Ch^\\u\\r+S,p,G + Ch^-^\\Ulu - Vhhp,G,h 

< Ch^\\u\\r+S,p,G + Ch^-''\\Ulu - u\\j,p,G,h 

+ Ch^~''\\u — Vh\\j,p,G,h 

< Ch^\\u\\r+S,p,G + Ch^-^\\u - Vhhp,G,h, (5.8) 

where IIJ) is defined as in the proof of Theorem 14.21 The inequality (15.81) means the 
following inequality holds 

Ch^-^\\u - Vh\\j,p,G,h > C||A'"^u||o.p,G - Ch^\\u\\r+S,p,G- (5.9) 

Together with Lemma l5.11 when h is small enough, we have 

Ch^-^\\u-Vhhp,G,h > C||A™^u||o.p.g = CA^IIuIIo.p.g > CX^. (5.10) 
In the second case, the integer i = 1 and r = m + 2m£. Then we have 

||V™A™^u||o.p.g = ||V™A™^(u - f;„)||o.p.G.h 

< ||V™A™^(u - n»||o.p,G.h + ||V™A™^(n> - f;,.)||o.p.G.h 

< Ch^\\u\\r+S,p,G + Ch^-^\\Ulu - Vhhp,G,h 

< Ch^\\u\\r+S,p,G + Ch^-^\\Ulu - u||,.p.G.h 

+Ch^~''\\u - Vh\\j,p,G,h 

< Ch^\\u\\r+S,p,G + Ch^-’"\\U - Vh\\j,p,G,h 

which means 


Ch^-^\\u-Vhhp,G,h > C\\V^A^^u\\o,p,G - Ch^\\u\\r+S,p,G- 
Similarly, together with Lemma l5.11 when h is small enough, we have 

Ch^-n|u-^;i*||,.p,G,/. > CIIV^A^^mIIo.p.g = CA^V^uHo.p.g 

> qG|A^||V™u||o.G > cx^+^ = CX^, (5.11) 

where the constant C depends on || V'"m||o g/|| V'"u||o o and the equality || V'"m||q n = 
Ai/2 

is used. 

Combining (I5.10p . (15.lip and the arbitrariness of Vh leads to 



vh&Vh 


'^h\\j,p,G,h 

h^-1 


> CX^. 


Together with the following property 

llu ||j,p,/i ^ II j,p,G,/i ^ ||f^ '^h\\j,p,G,h 
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we obtain the desired result (Oi) . The error bounds (EH) and EH can be proved 
similarly. □ 

Remark 5.1. The interior regularity result u G for a subdomain 

G CC n and S > 0 for problem a can he obtained from m Theorem 7.1.2] with 
f replaced by Xu. Now, we discuss some conforming and nonconforming elements 
that can produce the lower bound of the discretization error with the help of Theorem 
15.21 For the two-dimensional case (n = 2) there exist elements such as the Argyris 
and Hsieh-Clough-Tocher elements [5] that yield lower-bound results from Theorem 
15.21 for the biharmonic problem. The lower-bound results in this section hold for the 
Argyris element with m = 2, r = 6 (i = 1, £ = 1) and Ind6,used = {a ■ |q;| = 6}; and 
for the Hsieh-Clough-Tocher element with to = 2, r = 4 (i = 0, £ = 1), Indg^used = 
{a : |a| = 4}. 

Similarly, we can investigate the number of reliable numerical eigenvalues obtained 
from the finite element method for the biharmonic eigenvalue problem (in (|5.2I1 with 
TO = 2). The ingredients for the analysis are the lower bound result (15.911 and the 
asymptotic behavior for the j-th eigenvalue of the biharmonic operator. 

In [14], Pleijel proves the following asymptotic behavior for the n-th eigenvalue 
of the biharmonic operator: 



(5.12) 


For simplicity, we only state the following result and the proof can be given in the 
similar way of Theorem 14.41 

Theorem 5.3. Let N denote the total degrees of freedom in the finite element 
space Vh- Assume the mesh Th is quasi-uniform and we use the conforming finite 
element method to solve the biharmonic eigenvalue problem f id.21) with m = 2) and 
set je = for 6 G [0,1). Then the relative error of the jg-th numerical eigenvalue 
has the following lower bound 

and the absolute error of the jg-th numerical eigenvalue has the following lower bound 
Xj,,h - A,, > G(r - if<^rf^2rii-g)-4^ 


On the contrary, the number jg of the “reliable” numerical eigenvalues with relative 
error of Xjg converging at the rate h^O‘G~A) for 9 G (0,1], has the following estimate 

jg < C{r - (5 X5) 


Furthermore, the number jg of the “reliable” numerical eigenvalues with absolute error 
of Xjg converging at the rate g ^ following estimate 


jg < C{r - 


„{l_e){r-2) _ (l_e)(r-2) 

1)- ?- N -?- 


(5.16) 


6. Concluding remarks. In this paper, we have established some new upper 
and lower error bounds of the finite element approximations. As an important appli¬ 
cation, we derive lower bounds of the discretization error for the Laplacian eigenvalue 
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problem by piecewise polynomial approximation. We also show the asymptotic con¬ 
vergence behavior for the large numerical eigenvalue approximations. Our results 
reveals that the traditional upper bound is also a lower bound when we solve the 
eigenvalue problem by the finite element method. In particular, their dependence on 
the eigenvalue have the same power for both lower and upper bounds. This fact should 
have some illumination to the design of numerical methods for eigenvalue problems. 

Acknowledgments. The first author is supported in part by the National Nat¬ 
ural Science Foundation of China (NSFC 91330202, 11371026, 11001259, 11201501, 
2011CB309703) and the National Center for Mathematics and Interdisciplinary Sci¬ 
ence, CAS and the President Foundation of AMSS-CAS. The second author is sup¬ 
ported in part by the US National Science Foundation through grant DMS-1115530, 
DMS-1419040, and the National Natural Science Foundation of China (91430216, 
11471031). 


REFERENCES 

[1] R. A. Adams, Sobolev Spaces, Academic Press, New York, 1975. 

[2] T. Arbogast and Z. Chen, On the implementation of mixed methods and nonconforming meth¬ 

ods tor second-order elliptic problems, Math. Comp. 64 (1995), 943-972. 

[31 D. Arnold and G. Awanou, The serendipity family of finite elements, Found. Comput. Math., 
11 (2011), 337-344. 

[4] I. Babuska and J. E. Osborn, Finite element-Galerkin approximation of the eigenyalues and 

eigenvectors of selfadjoint problems, Math. Comp. 52(1989), 275-297. 

[5] I. Babuska and A. Miller, A feedback element method with a posteriori error estimation: Part 

I. The finite element method and some basic properties of the a posteriori error estimator, 
Comput. Methods Appl. Mech. Engrg. 61(1) (1987), 1-40. 

[6] I. Babuska and T. Strouboulis, The Finite Element Method and its Reliability, Clarendon Press, 

2001. 

[7] S. Brenner and L. Scott, The Mathematical Theory of Finite Element Methods, Springer-Verlag, 

2008. 

[8] P. G. Ciarlet, The finite Element Method for Elliptic Problem, North-Holland, Amsterdam, 

1978. 

[9] M. Crouzeix and P.-A. Raviart, Conforming and nonconforming finite element methods for 

solving the stationary Stokes equations. I. RAIRO Numer. Anal. 3 (1973), 33-75. 

[10] D. Gilbarg and N. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer- 

Verlag, Berlin Heidelberg, 2001. 

[11] P. Grisvard, Elliptic Problems in Nonsmooth Domains, Pitman Advanced Publishing Program, 

Boston, 1985. 

[12] M. Kffzek, H. Roos, and W. Chen, Two-sided bounds of the discretization error for finite 

elements, ESAIM: M2AN 45 (2011), 915-924. 

[13] Q. Lin H. Xie and J. Xu, Lower bounds of the discretization for piecewise polynomials, Math. 

Comp., 83(285) (2014), 1-13. 

[14] A. Pleijel, On the eigenvalues and eigenfunctions of elastic plates, Comm. Pure Appl. Math., 

3 (1950), 1-10. 

[15] R. Rannacher and S. Turek, Simple nonconforming quadrilateral Stokes element, Numer. Meth¬ 

ods PDEs. 8 (1992), 97-111. 

[16] G. Strang and G. Fix, An analysis of the finite element method, Prentice-Hall, Inc., Englewood 

Cliffs, N.J., 1973. 

[17] H. Weyl, Uber die asymptotische verteilung der Eigenwerte, Gott Nach. (1911), 110-117. 

[18] O. Widlund, Some results on best possible error bounds for finite element methods and ap¬ 

proximation with piecewise polynomial functions, Proc. Roy. Soc. Lond. A. 323 (1971), 
167-177. 

[19] O. Widlund, On best error bounds for approximation by piecewise polynomial functions, Nu¬ 

mer. Math. 27 (1977), 327-338. 

[20] Z. Zhang, How many numerical eigenvalues can we trust? http://arxiv.org/abs/1312.6773 

2013. 

[21] Z. Zhang, N. Yan and T. Sun, Superconvergence derivative recovery for the intermediate finite 

element family of the second type, IMA J. Numer. Anal., 21 (2001), 643-665. 


